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Polaronic effects are investigated in the two dimensional t-J model (describing Cu spins and holes 
of the planes in high-T c cuprates) in the presence of in-plane dynamic oxygen breathing modes. We 
focus on the effect of dynamic multi-phonon modes. Spin and charge structure factors are computed. 
We found evidence for self-localisation of the holes and increased short range antiferromagnetic spin- 
spin correlations between the remaining spins. This effect persists when all the in-plane oxygen 
vibration degrees of freedom are taken into account. This suggests that dynamic phonons could 
play a key role in stabilizing inhomogeneous stripe phases or domain wall structures. 
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The interplay of the electron-phonon interaction and 
localization, charge density wave (CDW) and spin den- 
sity wave (SDW) correlations in strongly interacting elec- 
tron systems is an interesting issue. In the static Born- 
Oppenheimer m — > oo limit, Zhong and Schiittler 
studied Hubbard and t-J models with an onsite Holstein 
electron-phonon apical oxygen coupling. They found on 
Hubbard (\/lO x ^/lO) and t-J (4 x 4) clusters that in 
the nearly half-filled case, in the presence of a strong 
onsite Coulomb repulsion and short range antiferromag- 
netic spin correlations, there were enhanced polaronic 
self-localization effects. In a recent paper A. Dobry 
et al. H investigated the effects of a single dynamic 
qo = (it, 7r) phonon mode coupled to a strongly corre- 
lated two-dimensionnal (2D) electron system modeled by 
a t-J Hamiltonian. For cuprates superconductors, this 
(it, 7t) mode corresponds to a planar breathing phonon 
involving correlated displacements of the in-plane oxy- 
gens. For vTO x y/lO clusters at one-quarter filling A. Do- 
bry et al. p| found evidence for a Charge Density Wave 
(CDW) instability at moderate values of the electron- 
phonon coupling. In addition, lattice deformations and 
CDW instabilities have been found in the one dimen- 
sional Holstein || and t-J-Holstein models. 

While both the Born-Oppenheimer and the dynamic 
single-mode approximations provide insight into this 
problem, one would like to examine the dynamic multi- 
phonon mode system. In particular, although the trunca- 
tion of the Hilbert space of the phonon mode is well con- 
trolled, the single mode approach of Ref. || nevertheless 
assumes that in-plane oxygens oscillate in phase. How- 
ever, there is experimental evidence that these modes 
have a small dispersion in the Brillouin zone, both in 
La.2- x Sr x CuOi or in bilayer Y BaiCu?,0§+ x materials 
In other words, the oxygen vibrational modes are 
very localized and can be considered as independent. In 
the following, we investigate the influence of taking into 
account a finite number M of oxygen modes by perform- 



ing exact diagonalisation studies of finite clusters. First, 
a simple model involving the coupling of a t-J Hamilto- 
nian with two independent qo breathing modes related to 
different sets of in-plane oxygen atoms is considered. Var- 
ious signatures of a tendency towards hole localization 
(which could lead to CDW or domain wall structures) 
and the enhancement of short-range antiferromagnetic 
correlations between the remaining spins are found. We 
reach similar conclusions on a 2 x 2 lattice when all the 
in-plane oxygen vibration degrees of freedom are taken 
into account. 

We start with a generalisation of the t-J-Holstein 
Hamiltonian JtJ, 
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where c\ a is the usual hole creation operator, n\ and 
are the electron and hole local densities respectively, 
m is the oxygen ion mass, fi is the phonon frequency 
and S = x, y differentiates the bonds along the x- and y- 
direction respectively. Throughout energies are measured 
in unit of the hopping integral t. Note that the electron- 
phonon g-term involves the coupling of each copper hole 
with the displacements of the four neighboring oxygens 
xi,s and Xi-s.s- For the purpose of our discussion it is 
convenient to re-write the electron-phonon coupling in 
momentum space, 
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where N is the number of sites, b A are the usual boson 
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creation operators, A q ^ 

The Hamiltonian of Ref. Q is obtained from Eq. (g) by 
(i) keeping only a single S, q = qo mode and (ii) replacing 
the sum ^= s by unity. This is a priori justified 
by the fact that the coupling constant A qi< 5 is maximum 
at wave vector qo. However, assuming in- phase oxygen 
displacements might overestimate cooperative effects of 
the phonons. On the other hand, a complete treatment 
of hamiltonian (2) requires to take the M=2N oxygen 
phonons into account. 

The first step to go beyond the single mode approx- 
imation (SMA) of Hamiltonian (j^) is to consider two 
independent bosons & qo , x and & qo , y (M=2) for the hor- 
izontal (5 = x) and vertical bonds (5 — y). In this 
case, the sum s is replaced by ^ ^ 5 and only 

two q = qo modes are retained. The ground state (GS) 
energies Eq(Xq) of the single and double mode approxi- 
mations (DMA) of Hamiltonian (Q) are calculated on a 
N = \/20 x V20 cluster doped with two holes by trun- 
cating the (strictly speaking infinite) Hubert space of the 
phonon mode(s) to a finite number of bosonic levels N p h- 
This can be achieved, for example, by keeping only a 
finite number of levels n p h + 1 for each oscillator. In 
other words, a maximum number of n p h local phonons 
could be excited for each oxygen mode. Then the to- 
tal number of degrees of freedom for the bosons becomes 
Nph = (n p h + 1) M where M is the number of modes re- 
tained (M=l for SMA, M=2 for DMA and M=2N for the 
full problem). 

The GS energies E (\ )-E (0) for M=l and M=2 are 
displayed in Figs. |l](a,b) in a regime (J = 0.3) where, 
in the absence of phonons, spin fluctuations leads to a 
small effective attraction between the holes [||. The de- 
pendence on the number of levels n p h + 1 is shown in 
Fig. [1(a). For M=l, as in the quarter filled case con- 
sidered in Ref. [||, the truncating procedure is well con- 
trolled and the GS energy decreases with N p h and typi- 
cally converges for N p h ~ 20 or 30. When more than a 
single mode is included, it is clear from Fig. |l|(a) that 
n p h (and not N p h) is the key parameter controling the 
convergence of the GS energy. We then expect that a 
similar number of levels per phonon mode (~ 15) is still 
required to fully account for the vibrational degrees of 
freedom when M > 1. Hence, disk space limitations 
clearly restrict the study of the full hamiltonian (2) (with 
M = 2N) to small 2x2 clusters (see below). However, 
useful informations can still be obtained with n p h ~ 2- 
5 since most physical properties are smooth monotonic 
functions oin p h- Fig. |l|(a) suggests that polaronic effects 
are stronger when the number of independent modes (M) 
is increased. As seen from Fig. |l|(b), for increasing Ao 
both models show a crossover from an adiabatic regime 
(the energy depends weakly on the coupling strength) to 
a polaronic regime. We, indeed, expect that, for a suffi- 
ciently large value of the electron-phonon coupling, the 



system will start to form large hole polarons. By com- 
paring the set of data for the two models we observe that 
the onset A co corresponding to polaron formation might 
not depend too strongly on the truncation parameter M 
(or n p h) although for Ao > A co polaronic effects increase 
with M. 

In addition to the formation of polarons, the electron- 
phonon coupling (treated in the SMA) also favors the 
configurations where the holes (polarons) stay dynam- 
ically on the same sublattice. Such a polaron-polaron 
interaction can eventually lead to a CDW instability as, 
for example, it is found at quarter filling [g] in the SMA 
or to domain walls formation. To complete our study, we 
address now the issue of the spin and charge correlations 
at low hole doping (10%) and study the influence of the 
number of phonon modes on the (tt, n) structure factors. 
We shall consider here only the case of commensurate 
(tt, tt) charge fluctuations since the density is too close to 
half filling and the system probably not large enough to 
accomodate any incommensurate order. The charge and 
spin structure factors at momemtum (tt, tt) are given by 
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where |^o) is the GS wavefunction calculated both in the 
SMA and DMA cases. 5 c (qo) is shown in Fig. |2| (a) as 
a function of the electron-phonon coupling. For a single 
mode (SMA) a significant increase of the charge structure 
factor occurs around Ao = 0.45 indicating a tendancy to- 
wards CDW formation or at least towards some localisa- 
tion of the charge carriers. Interestingly enough, Fig. || 
(b) shows that this transition gives rise to a simultaneous 
increase of the spin-spin antiferromagnetic correlations. 
This phenomenon might be related to a sudden increase 
of the polaronic mass or to the localisation of the holes in 
some kind of (fluctuating) CDW state or stripe (domain 
wall) phase. Less mobile holes clearly lead to a weaker 
perturbation of the antiferromagnetic background, hence 
leading to enhanced short range spin-spin antiferromag- 
netic correlations. It should be noticed that, on the con- 
trary, for quarter filling Q the CDW transition is fol- 
lowed by a decrease of the spin correlations. When two 
oxygen vibrational modes are considered (M=2), Fig. ^ 
(a) and (b) show that the increases of the charge and spin 
structure factors (indicating some hole localization) oc- 
cur for approximately the same couplings than for M=l. 

We finish our study by the full treatment of all the in- 
plane oxygen degrees of freedom as seen in Fig. |^(a,b). 
Although the system size is smaller in this case (2 x 2) 
the results are qualitatively similar. Fig. ||(b) clearly 
shows a large increase of the spin structure factor at A co . 
This increase becomes sharper and sharper when more 
phononic levels arc included (eventhough the GS energy 
is not fully converged, see Fig. ||(a)). Note that the cross- 
over value seems to be smaller for n=0.75 (A co ~ 0.3) 
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than for n=0.9 (A co - 0.4). 

The previous results can be physically interpreted as 
follows; In the SMA the oxygen atoms are constrained to 
oscillate in phase leading to a dynamical potential which 
favors the occupation by the carriers of one of the two 
sublattices. As discussed in Ref. (|] this can lead to a qo 
CDW instability as can be seen for example in a mean- 
field treatment. The DMA is the simplest way to intro- 
duce some dephasing between the harmonic oscillators. 
In the extreme case of out-of-phase motions of the two 
oxygen sets the effective potential for the holes exactly 
vanishes. Therefore, one would naively expect that the 
DMA or the complete treatment of hamiltonian (2) lead 
to a weaker tendency towards CDW than the SMA. On 
the contrary, our results show that, in these more reliable 
treatments, the electron-phonon coupling still leads to 
strong self-localization effects || . More importantly, they 
also suggest that dynamic in-planephonons could help to 
stabilize inhomogeneous phases where holes are 

squeezed into stripes or clumps and where the remaining 
antiferromagnetic regions have enhanced S s (ir, if) corre- 
lations. 
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FIG. 1. Pair energy calculated on a v20 x \/20 t-J cluster 
with 2 holes for J = 0.3 and Q, = 0.2. Open (solid) symbols 
correspond to a single (double) (it, it) breathing mode, (a) 
Energy vs number of phononic levels included in the calcu- 
lations, (b) Energy vs coupling strength Ao. The number of 
phononic levels considered is indicated on the plot for each 
case (single and double modes). 

FIG. 2. Charge (a) and spin (b) structure factors at mo- 
mentum (it, 7t) vs electron-phonon coupling calculated on a 
\/20 x -^20 t-J cluster doped with 2 holes for J = 0.3 and 
fl = 0.2. Open (solid) symbols correspond to a single (dou- 
ble) (n, 7r) breathing mode. N P h = 25 phononic levels were 
kept in the calculation. 

FIG. 3. GS energy vs n p h (a) and spin structure factor at 
momentum (n, n) vs electron-phonon coupling (b)calculated 
on a 2 x 2 t-J cluster doped with 1 hole for J — 0.3 and 
Q — 0.2. The total 8 independant oxygen modes have been 
considered. 
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